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I, INTRODUCTION 


The Material Balance Method of calculating the behavior of pet 
roleum reservoirs is one of the petroleum engineer's most important 
tools. Introduced in 1936, the Material Balance has undergone many 
improvements intended to simplify the mathematical techniques of handling 
the complex physical relationships involved. In spite of all advances 
in technique, the method remains unwieldy and laborious. 

The basic concept of the Material Balance applies the Law of 
Conservation of Mass to the process of producing petroleum from an under- 
ground reservoir. The amount of hydrocarbons originally present in the 
reservoir must equal the amount produced plus the amount remaining in the 
cround at any given time during the production life of the reservoir. As 
normally used, the Material Balance is written volumetrically. It relates 
the volumes of oil and gas (produced and remaining) by means of their 
physical characteristics (gas volume, gas solubility and oil shrinkage) 
at the pressure prevailing in the reservoir. 

The Material Balance has been most widely used for the solution of 
two general problems: the estimation of original oil in place and the 
prediction of future production performance of a reservoir. The reservoir 
discussed in this paper will be of the Solution Gas Drive Type, with fixed 
reservoir volume, and having no gas cap and no water encroachment into 


the reservoir. 


A. Estimation of the Original Oil in Place 


The quantity of oil and gas originally present in a reservoir may 
be estimated by observing the drop in reservoir pressure during the pro- 
duction of a certain amount of fluids, and relating these data to the 
laboratory-measured behavior of the reservoir fluids caused by pressure 
changes. The volume of reservoir fluids (with known shrinkage and 
solubility characteristics) which will undergo a certain pressure change 
because of the removal of a given quantity of these fluids is a fixed 
quantity. Thus, the drop in reservoir pressure caused by a given amount 
of oil and gas being produced defines the volume of these fluids origi- 


nally present. 


B. Prediction of the Future Production 


Performance of a Reservoir 


For any given pressure assumed to prevail in the reservoir at some 
future time, the volumes of oil and gas produced (when related to their 
physical characteristics at that pressure) constitute a unique solution 
to the equation balancing the original volume of fluid with the volumes 
produced and remaining. An additional concept is involved however. 

The permeability of the porous reservoir medium to oil and to gas changes 
as production proceeds. Thus, as the saturation of liquid in the reser~ 
voir decreases as oil is produced, the relative amounts of oil and gas 
which flow to the well will change. This phenomenon cannot be directly 
related to the Material Balance itself, yet it determines the relative 
quantities of oil and gas which will be produced. The relationship be- 
tween the relative permeability to oil and gas and the liquid saturation 
in the reservoir must be known or assumed, since it must be considered in 
arriving at the amounts of oil and gas vc be produced. These amounts of 
gas can then be used to predict a solution to the Material Balance 
at some assumed future point of pressure and production. 

It will be the purpose of this thesis to develop a graphical means 
for predicting the performance of a depletion drive reservoir based ona 


modification of the Material Balance equation. 


II. REVIEW OF THE LITERATURE 


The use of the Material Balance for estimating oil in place was 
developed by Schilthuis in 1936.7 His equation (shown here for a solution 
gas drive reservoir without gas cap, water encroachment or water production) 


niu (ey, = 8, )¥ | . 


occa (1) 
still remains the most used basic form of the Material Balance. 

At about the same time Katz- proposed a tabular method for evaluat-— 
ing oil in place. This method was later shown by Pirson> to be equivalent 
to the Schilthuis method. 

The use of the Material Balance Equation for orediction of future 
reservoir performance was proposed considerably later. Babsan" developed 
a trial and error solution based on the Material Balance Equation and the 
Instantaneous Gas-Oil-Ratio Equation. His method was very cumbersome and 
required a great deal of computation. It is rarely used tcday. In the same 
year, 19h, pseuee proposed a solution to the problem which was consider- 
ably simpler to use. : 

Nuskat® developed a form of the Material Balance expressed as a dif- 
ferential equation and applied it to prediction of depletion drive reser- 
voir perfornance. 

The most widely used of the Material Balance prediction techniques 
is the method using a trial and error solution of Schilthuig’ Equation!» 8 


This method requires the simultaneous satisfaction of three equations: 


References are listed in the Bibliography. 
Nomenclature shovwm in Apper.dix I. 
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The Material Balance Equation 


n{u + (e ~ s)v] 


= ee (1) 
U =U, 
The Instantaneous Gas-Oil-Ratio Equation 
Mo BK 
rT: = No B Ke + § (2) 
L Plow K. 


The Liquid Saturation Equation 


Sy (lee Ca (3) 


The trial and error solution may be conducted in the following 


manner s 


Assume a certain drop in reservoir pressure 

Estimate the volume of oil which would be produced during 

this pressure drop 

Using this estimated value of oil production, solve the 

Liquid Saturation Equation, (3), and obtain the corresponding 
K 

value of i from the known or assumed ze vs. Sy relation= 

Ship. : 

Solve Equation (2) for the value of Instantaneous Gas-Oil- 

Ratio which would obtain at the assumed pressure. 

Using the assumed and estimated values of 011 production, gas- 

oil-ratio and fluid characteristic solve the Material Balance 

Equation. If a balance is obtained, the assumed values were 

correct. If the equation does not balance, a new value of oil 


production or pressure is estimated and another trial is made. 


This process is repeated until the equations are satisfied at sufficient 
pressure points to define the future performance of the reservoir. 

The use of the Schilthuis Equation, while more flexible in the 
handling of such conditions as water encroachment and gas cap expansion, 
is extremely tedious. For the simplest reservoirs, twenty~five columns 
of calculations must be used, and five or six trials at each of seven or 
eight pressure points is not unusual. All this must be performed with a 
calculating machine. Despite this, it is still the most widely used 
mathematical prediction technique. 

In January, 1955, a prediction method based on the Schilthuis 
method, but greatly simplified, was published by Tracy.” The Tracy method 
is the single biggest advance in the use of the Material Balance Equation. 
This method will be described for the case of a depletion drive reservoir 
without gas cap or water encroachment. 


Tracy took the Schilthuis Equation 


nlu ar, = 3y)v] 


i (1) 
and substituted values for u, u, and ry, 3s follows: 
u= 8 + (85-8) v (i) 
uy = 8, (5) 
mmo S (6) 


This substitution resulted in; 


ES 


(e 


(7) 


aye 


which may be expanded to: 


(74 \ 
" niy- s} He G 
N ja Zo \ (ee _ = = fo . _ \ (8) 
vo s ) igo. Kv / v o/ 


Thus, n and G are multiplied by coefficients which are functions of pressure 


only. These coefficients are: 


¢ - s) 
‘ ‘Ere a) i 


ps, 2) 
é ) 
= (10 
© & -2) eee 
Vv y oo G} 
Equation (8) then reduces to: 
Nos ng, + 64 (12) 


In this form, the equation may be used to estimate original oil in place. 
Tracy's prediction method uses this equation in two forms with N 
taken as equal to one barrel of stock and oil. 


1. As the Prediction Equation: 


1 = (nya 6, + Gyn 6g) 
An = (12) 


Ns +7 
a i=l 
a 


2. As the Material Balance Equation: 


Orhan + 8p bas 1 (13) 


Equations (12) and (13) are used in conjunction with Gas-Oil-Ratio 


Equation: 


ol oe (2) 


and Liquid Saturation Equation: 


Sms a9 (Pos (3) 
Bo 


Equations (2) and (3) being linked by knowledge of the a vs. Sy; relation- 
ship. 

The use of Tracy's Prediction Method calls for the following 
procedure 

1. Starting with a point in production history (with known values 
of ns, G, 6; re for a given pressure), assume a new lower pressure point. 

2. Estimate an instantaneous gas~oil ratio for this pressure 

3. Using this pressure and gas-oil ratio in the Prediction Equation 
(12), compute a value for An. 

4, Add this value of Mn ton, -1 (n, from the previous pressure) 
to get the value of n, for the new assumed pressure. 

5. Solve the Liquid Saturation Equation, (3), using this value of 


n, and obtain a value for 


° 


ober 


6. Using the ze value, solve the Gas-Oil Ratio Equation. This 
fe) 
value of Instantaneous Gas-Oil Ratio should equal that estimated in step 2. 
If it does not, use the value of cy found in step 6 as the estimated value 


in step 2 and re-do steps 2 through 6. When the values of Instantaneous 


Gas-Oil Ratio in steps 2 and 6 are equal, substitute tne values of ns and 


Gs into the Material Balance Equation (13) and solve. The sum of 
ng Py + G, 4, should be 1.00 + 0.02. 

The number of trials required is cut down to two or three in Tracy's 
method because Instantaneous Gas-Oil Ratio - a relatively insensitive 
factor - is estimated. This cannot be done in the Schilthuis Equation 
as the cumulative value of gas-oil ratio, rather than the instantaneous 
value, is used in the equation. The cumlative gas-oil ratio value cannot 
be estimated with any degree of accuracy. 

Sturdivant, 10.11.12 in a series of three articles, expanded upon the 
use of the Tracy method developing a formula to be used in the event of gas 


reinjection into the reservoir. It was also shown that accuracy to three 


significant figures in the calculation was sufficient. 


IIIT. STATEMENT OF THE PROBLEM 


This investigation was conducted with the object of developing a 
graphical method for solution of the Performance Prediction problem. It was 
felt that the Prediction Equation developed by Tracy, with its separation of 
the pressure factors into coefficients, would lend itself to such treatment. 
In the Schilthuis Equation, the complexity of form appeared to preclude any 
simple graphical treatment. 

Of the various graphical problem solving methods available, nomography 
appeared as the best for trial and error soluticns. The forms of the equa= 
tions indicated that nomographs of reasonable simplicity could be constructed 
to solve them, For trial and error solutions, it was felt that the nomo-— 
graphs should be simple enough that the operator could see at a glance how 
much a given change in one factor would affect the other factors in the prob= 
lem. Also, with simple, one-setting nomographs, the operator can readily 
enter the diagram with the dependent variable (or answer) and work backwards 
to find the value of indevendent variable (or normal input). This is of value 
on trial and error solutions. For these reasons, as well as for ease of con- 
struction, it was desirable that the nomographs be of the type in which a 
single line crosses all scales, obtaining the answer in a single setting. 

There are a number of nomographs in the literature concerned with 


13514515516 All of them are 


various aspects of the Material Balance problem 
quite complex, requiring entry with several individual pressure variables, 
and requiring several settings of a line before a solution is obtained. 


These nomographs are general in application and could be used on reservoirs 


having the same ranges of variables as the nomographs. 


10 


ales 


By constructing nomographs for a specific reservoir, ali of the fluid 
characteristics (which are functions of reservoir pressure) and combinations 
of them, can be treated as a single, or recurring, variable of pressure. 
This reduces the mimber of variables to the point where simple, single] 
setting nomographs can be constructed. Construction of the nomographs for 
a single reservoir allows selection of optimum scale lengths to specially 
handle the ranges of the variables for that reservoir. More generalized 
nomographs would require the ranges of the scales to be broad enough to fit 
all or several reservoirs. This would result in less accuracy for any 


particular reservoir. 


IV. METHOD OF INVESTIGATION 


The four equations used in the Tracy method of prediction were 
separately investigated to find the forms most likely to yield the simplest 
nomographs. These equations, given previously as Equations (12), (13), (2) 


and (3) will be considered separately. 


eat 
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A. The Prediction Equation 


: ; Pest wt 
This equation was given ? ase 


An = a (14) 
é, +7, 8, 
where 
re + Ye 
ee raed (15) 


It is readily rearranged into this form: 


(n,, + 4n) é. +16, a * te Zon) b, =a (16) 


In this form, it is seen to be equivalent to the Material Balance 
Equation as; 
ns, + An = ns (17) 
and 


Gaky te, 2 AG (18) 
If these values are substituted into Equation (16), the result is the 
Material Balance Equations 


n, 6 +G. g =1 (13) 
1 sn eS 


Returning to the rearranged form of the Prediction Equation: 


(n. 


oe Lyn) g + (G, +r An) do =1 (16) 


=] a “g 


it is noted that there are four variables: 


(n, 4 + fn), a function of oil production 
Let this be a (u} 
(Ga + are ee n), a function of the gas produced with the above 
amount of oil 
Let this be f, (v) 
ae a function of reservoir pressure 
Let this be f, (Ww) 
B.,3 a different function of reservoir pressure 
Let this be f), CW) 


Thus, Equation (16) can be written as: 
fy (u) x f, (W) + £, (v)s fy, (WW) = 1 (ry) 


Note that the pressure, represented here by "W" is a recurrent 


variable. 


Divide Equation (17) by f, (W), then: 


i (u) ot £5 (v) 2 F W = #> (Wy (18) 
| 3 

Now let 
fy CW) 
ee S| | LT (N) a different function of pressure 
f3 (W) 5 

and 
1 oe 

FW) = f, (W), another function of pressure. 

Equation (18) now becomes: 
f, (u) + £5 (w)rty (W) = £, (W) (19) 


a 


This form is capable of representation by a simple nomograph having 


the following configuration: 


n scale G scale 
(ny = D447 La) (G. = G4 +On ee) 
| id 
Enter Read | /Enter 


Tne geometric proof for this nomograph is given in Appendix IV. 


Now, evaluating the pressure functions fo (W) and fi (W) results ins 


f, (WwW é. 
ie oa ees Looe. (20) 
af 3 (w} é [f 7 a 
rg “) 
(@o 


al 
W a ig : fat me ws 
f¢ (W) on, - aoe 


2 


These values, when substituted into the Prediction Equation, gives 


(G. 4 uae a n) (Po ~ § \ 
: e SF xe ta rw) 


f pe Oy Re a Ee a 
(n, . + An) ? : z (22) 
i-l eo | a & 
vv F , 7 ow Ss 


Equation (22) is the form of Tracy's Prediction Equation which is 


solved by the Prediction Equation Nomograph,. 


te we 
- 7 ’ 
‘en aeee com eee eT eet te 
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B, The Material Balance Equation 


This equation is equivalent to the Prediction Equation and can use 
the same nomograph. One solution of the nomozgraph serves to satisfy both 
the Prediction Equation and the Material Balance Equation. 

Placed in the form used to construct the nomograph, the Material 


Balance Equation becomes: 


eae cm eo) 
OS i en ae sa 
Vv Lv 


Equation (23) is solved by the Prediction Equation nomograph. 


PaT ei jee sist 
ai 4 
S'* (ae seaeeteS aS 


4 in? wad [he 


yi 
C. The Liquid Saturation Equation 


The Liquid Saturation Equation is used to obtain a value of =® from 


ee K 
a known and extrapolated relationship between Sy and _8, 


(8) 


In the equation 


5.2 Se oes Se (t-\g 


The equation can then be rearranged: 


Sy -K, = (K,) (8) (i - ay) (2h) 


Log (Sy; + K,) = Log Ky + Log (1-n,) (25) 


Equation (25) was used to construct a parallel line nomograph with 
logarithmic scales having the following form: 


Log (l-n,) Log (S,~K,) Log (Ky B) 


on ee oe 


*For proof of the geometry of this nomograph, see Appendix IV. 
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Log (1-n, ) is a simple function of n, and the left scale was therefore 
graduated directly in units of n,. Similarly, ff , and therefore Log Ky 
is a function of pressure, so that the right scale was graduated in units of 
pressure. The center scale was graduated on one side in units of Sr, and by 


use of the known ze vs. S; relationship, the other side was graduated in 


fo) 
units of ze - The finished nomogram appeared as follows: 
) 
P Fe 
bt, 
Fle ey 3 | 
7 ‘Enter 


19 
D. The Gas-Oil Ratio Equation 


In the instantaneous Gas-Oil Ratio Equation 


Mo O@ ky, s (2) 


iy ane ee 


Mo a » V, and s are specific functions of the reservoir pressure. 
By combining all of these but s into a single function of pressure, the 


equation was rearranged to: 


K 
ae a Roth ue s8 (26) 


This form is identical to 


sii (u) ~ f,.(v) ° ts (W) = f, (W) (27) 


h 
where = (W) and fi, (W) are different functions of the recurrent variable, 


reservoir pressure. 
A nomograph, similar in theory” to the Prediction Equation nomograph 


was constructed. The finished form appeared as follows: 


My 
“e p |% 


*See Appendix IV for geometric proof, 


aT 
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E. Average Gas-Oil Ratio Equation 


It was found desirable to construct a very simple nomograph to solve 


the equation 


re Shae 
eS ya (15) 


The equation was rearranged to: 


A nomograph with uniform parallel scales set equal distances apart 


was constructed. 


Entry can be made on any two scales and a value read on the third 


scale, 


V. USE OF THE NOMOGRAPHS TO PREDICT THE PERFORMANCE 


OF A SPECIFIC RESERVOIR 


To test the application of the nomographic prediction method, an 
actual reservoir was selected and nomographs constructed using the data 
available concerning this reservoir. The particular reservoir chosen was 
a mid-continent limestone reservoir, which is still in the early stages of 
production. For this reason, it will be called Reservoir "X." The basic 
data available on Reservoir 'X" are listed in Appendix V. In addition, two 
sets of prediction curves (Pressure and Gas-Oil Ratio Vs. Cumulative 
Production) were available for comparison with the results obtained in this 
investigation by the nomographic method. One of these sets of prediction 
curves was obtained by trial and error use of the Schilthuis Equation. The 
other was obtained by using Tracy's method. 

The insensitivity of instantaneous gas-oil ratio values over a range 
of incremental oil production values was mentioned in connection with the 
Tracy method. The relative positions of the n, G and P scales on the 
Prediction Nomograph show this feature graphically. The pressure scale is 
very close to the gas production scale. For this reason, it is best to enter 
the nomogram with an assumed value of oil production, rather than to assume 
a gas-oil ratio as in the tabulated Tracy method. 

The following approach was found to give the easiest solution using 
the nomograms: 


1. Set up a tabulation sheet as follows: 


ed 
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(0) P 1209 1100 900 
(A known point of 
Production History) 


(1) An 0.00773 0.0128 
(2 Seo 020319 0.03963 
(3) ny 0.0319 0.03963 0.05211 
(4) Glen aubr: 46.75 65.5 108.3 
(5) Ginn 46.75 _ 65.5 
(6) rg An 18275 42.8 
(7) ra 21,25 34130 
(8) Ts] 21,00 24.50 
(9) rj 21,00 2450 420 
(10) Kg 0.0362 0.085 
(11) oe (From G.0.R. 21,50 11,20 
+ Equation) 

(12) (nj) (1.8 x 10") Toto 9375 
(13) +575 ,000 575.0 575.0 
(1L) Cumlative Production 1288.0 1512.5 


(in 103 bbls.) 


2. Select a point of actual production history, knowing pressure, 
oil production from the bubble point to that pressure, n;, expressed as a 
fraction of original oil in place, gas production from the bubble point to 
that pressure expressed as a function of original oil in place, and instan- 
taneous gas~-oil ratio, Tye 

3. Assume a new pressure, 100-200 psi lower than the historical 
point used. The nj, G; and r; of the historical pressure point now become 


Ns js Gs_4 and P54 in the above tabulation, 


lh. Estimate an increment of oil production, An. Enter this on 


23 


line one of the tabulation. Add to it the nea The sum is n,- 

5. Using this value of n; at the assumed pressure, erter the pre- 
diction nomograph on the n and P scales and read the value of 
Gij+ Anr, onthe G scale. Enter this value on line ) of the 
tabular sheet and subtract the Go from it. The difference is An re 

6. With a slide rule, divide the value of An Page (line 6) by 
An (line 1) to get r, (line 7). 

7. Enter the average gas-oil ratio nomograph with r;_, and e and 
find rj (line 9). 

8. Using ny (line 3) and Pressure, find Zz from the liquid 
nomograph and enter this value on line 10. ° 

9. Enter the gas-oil ratio nomograph with ze and Pressure and 
find the r; value based on the gas-oil ratio ae rane Place on line ll. 

10. Compare lines 9 and 11 which will be equal if the assumed value 
of An was correct. 

Note that as assumed values of Z\n are increased, resulting values 
of r; on line 9 (from prediction nomograph) decrease sharply, whereas values 
of yr, 6m line 11 (from gas-oil ratio nomograph increase slowly. Thus, if the 
r; of line 11 is less than that of line 9, a slight increase of An will 
bring them together. If the two values of r; are within a few hundred of 
each other, set the An r, value from line 6 on the slide rule "D" scale, 
place the value of = (line 7) opposite it on the "C" scale, and read a 
nearly-correct value for An under the index on the "D" scale, 

If this new value of An causes a radical change in n,;, it may be 
necessary to enter the nomograph for a new Gs4 + 25 1 roe However, from 
the scale positions, it can be seen that Gii4 ae ae a De is very insensitive 


to changes in An at most pressures. 


2h 


11. If the new value of Mn caused more than a slight change in 
ni» re-enter the liquid saturation and gas-oil ratio nomograms to obtain a 
corrected rs; for line ll. 

12. Lines 9 and 11 should now be within perhaps 50 or 60 of each 
other. Alter An slightly to bring line 9 in exact agreement with line 
10. Add the corrected On ton;_,. The result is n,. 

13. Multiply n; by the volume of original oil in place for line 
12 and add line 13, the production from the original pressure to the bubble 
point. Line 1h is the predicted cumulative production at the assumed 
pressure. 

1). Assume a new pressure about 200 psi lower and repeat the 
process. 

Special care in the averaging of r; and rs_; must be taken when 
computing at about the peak of the Gas-Oil Ratio Vs. Production Curve, as 
any averaging process assumes local linearity of a curve. It is best to 
use Prone wiasears increments in this region. 

Using the method described above, about two hours were required 
to calculate the predicted performance using seven pressure steps. 

An alternate method was developed which takes about the same time 
to use, but which may be a little more accurate: 

1. Use the same tabulation sheet as shown on page 22 above. In 
addition, a large "scratch graph sheet" of squared paper (10 x 10 to the 
inch was used) and a sheet of scratch paper are necessary. 

2. Same as step 2 above. 

3. Same as step 3 above. 

h. Estimate three values of An quite close together. Jot on 


scratch sheet or carry mentally. 
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5. Enter the prediction nomogram with the middle value of An 
and read a value of Ge + Anr,- (It will be seen that the change in 
Gaede yan r, over the three values of Mn is so small as to be un~ 


readable.) 


6. Subtract G,_, from G,,+ Onr,. Divide An r, by each of 


-1 
the three values of An from step obtaining three values of rae 

7. On the scratch graph sheet, plot the three values of r, VSe 
the respective “An values. 

8. Add the three values of An to n,_,- Use the resulting 3 values 
of n; in the liquid saturation anc gas-oil ratio nomographs to obtain three 
values of rj. 


9. Using the average gas-oil ratio nomograph and entering with 


eas 


vo and the three rs values of step 8, find three values of roe 


10. Plot these values against An values on the scratch graph 
sneet. 

11. Where the lines of step 7 and step 10 cross, read the correct 
values of gvnand rag. Enter these values on the tabular sheet and check 
the computation on the nomographs, comparing tabulation lines nine and 1]. 

12. Make minor adjustments of Mn to match line 9 to line ll. 

13. Same as steo 13 above. 

1h. Same as step 1) above. 

Figure 1 is a graph comparing the results obtained by the nomo- 
graphic method with those obtained by trial and error solution of the 
Schilthuis Equation and by Tracy's method. The results in general are 
in good agreement. The higher gas-oil ratio values shown on the Schilthuis 
method curve are thought to have resulted from use of different values of 
"vy" than were used in the Tracy method and in the construction of the 
nomographs. The Schilthuis method curve was taken from an engineering 


consultant's report on Reservoir "X", 
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VI. SUMMARY AND CONCLUSIONS 


This investigation had as its object the development of a graphical 
technique for predicting production performance of a petroleum reservoir, A 
graphical technique involving the use of nomographs based on the Tracy method 
was worked out. Using this method nomographs were constructed to predict the 
performance of a specific actual reservoir called Reservoir "X." Tha per- 
formance of Reservoir "X" under depletion drive was predicted using these 
nomographs., The results, compared with the results obtained by other 
methods, show that the use of nomographs constructed for a specific reser~ 
voir gives a forecast of reasonable accuracy. 

The total time involved in nomograph construction and in solving the 
prediction problem with these nomographs proved to be about the same as that 
required to solve a reservoir problem by the Tracy method. However, reser= 
voir performance predictions are often repeated, starting with a later point 
in production history, as the field develops. With this in mind, it is be- 
lieved that the nomographic method would save considerable calculation, A 
repeated prediction could be run in about two hours using the same nomo- 
graphs, A repeated calculation, even by the Tracy method, would take #:-% 
considerably longer. 

As the production history of the field develops, it may be found that 
the basic field data relating relative permeability ratio and saturation were 
extrapolated in error, In this event, new and more correct values of ze 
may be marked on the center scale of the liquid saturation nomograph 
without affecting its accuracy for subsequent use, 

The limited time available precluded investigating the application 


of this technique to a reservoir undergoing gas reinjection. It is almost 
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certain that a set of reasonably simple nomographs could be constructed to 
handle predictions under various conditions of gas reinjection. 

The time available was also too short to allow the method to be ap= 
plied to reservoirs of differing characteristics. In particular, it should 
be tried on a sandstone reservoir, on a reservoir with a high initial pres- 


sure, and a reservoir with a low initial pressure. 


1, 


hy 
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APPENDIX TI 


Nomenclaturs 


Me Vv 
Cumulative gas production, standard cubic feet 
G at reservoir pressure being considered 


G at previously considered reservoir pressure 


Relative permeability ratie, a dimensionless number 


otal original oil in place, in stock tank barrels 

Cumulative oil production, expressed either in stock tank barrels 
(or as a fraction, when N is taken as equal to one barrel) 

The value of n at the reservoir pressure being considered 

The value of n at the last previously considered reservoir 
pressure 

Average gas~-oil ratio, standard cubic feet per barrel of stock 
tank oil 

Instantaneous gas-oil ratio at reservoir pressure being con~ 
sidered standard cubic feet per barrel of stock tank oil 
Instantaneous gas-—oil ratio at last previous reservoir pressure, 
standard cubic feet per barrel of stock tank oil 

Cumulative gas~-oil ratio, or the ratio of total standard cubic 
feet of gas produced to total barrels of stock tank oil 
produced 

Total liquid saturation in the reservoir 

Oil saturation in the reservoir 


Saturation of connate water in the reservoir 
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Solution gas-oil ratio, or the solubility of gas in crude oil 
at reservoir pressure and temperature, cubic feet per barrel 
Value of s at original reservoir pressure 

Volume in the reservoir occupied by one barrel of stock tank 
oil plus all the gas originally in solution in that oil, or the 
two-phase formation volume factor 

Value of u at original reservoir pressure Us =6 6 
Barrels of free gas space occupied in the reservoir by one 
standard cubic foot of gas 

Formation volume factor, or the volume occupied at reservoir 
conditions by one barrel of stock tank oil plus its dissolved 
gas 


Value of # at original reservoir pressure 


Increment of oil production between two reservoir pressure 
steps 

Viscosity of reservoir oil at reservoir conditions, centipoises 
Viscosity of produced gas, at reservoir pressure and temperature 


centipoises 
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APPENDIX III 


Considerations Affecting the Accuracy of Nomographs 


1. Dimensions of the Diagram 


The accuracy of a nomograph can be increased by increasing the 
dimensions of the diagram, within reasonable limits. This is due to the 
increase of relative accuracy and fineness of scale graduations. 

2. Order of Accuracy 

Nomographs, in general, have a degree of accuracy exceeding that 
of a slide rule of similar dimensions. In the nomograph, the ranges of 
the scales are tailored to the ranges of the variables involved. Short 
ranges are often expanded into long (and, therefore, more accurate) scales. 
3. Decrease in Accuracy due to Complexity 

The accuracy of a nomograph decreases sharply when more than one 
setting of a line must be made to obtain the answer. The more inter- 
related lines that must be drawn on a given nomograph to obtain a solution, 
the less accurate the answer will be. 
lh. Relative Placement of Scales 

The accuracy of a nomograph is affected by the relative positions 
of the scales. If possible, the scale on which the answer is to be read 


should be between the other two scales. 


APPENDIX IV 


Geometric Proof for Nomographic Forms Used 
A. Prediction Equation and Material Balance Equation 
The equation form for the nomograph is: 


f, (a) et or) f. (NW) = uy CH) (19) 


2 


Take two parallel axes for the variables u and v, anda 


curved axis in general for the variable WV: 


The line AB, of length k, joins the zero points of the scales for u and v. 
Points x and y are any points on the u and v scale, x = my fy (u) and 
y = mp fo (v), where mM, and mp) are the respective scale moduli. 

Take any index line CD intersecting the scale for W at BE. Con- 
struct EF parallel to the u and v axes. Let EF = z and AF = z'. From 


Similar triangles SCH and EDG: 


ou 


X-Z 2 2-y 
vel k-z! 
(k-z')x = kz -— 2z' + gz! — yz! 


(k-2')x + ya! = kz 


y 7% ‘ eee _ ei k ee 
care ies 
Substitute for x and y: 
My Pe) oF a \{ f, (v)j= ae 
teal en esse oe ieee ee 


This equation will become the original equation if fy (W) and f¢ (W) have 


the following values: 


My% 
m, (k-2") 


fy (W) 


k 


fe (W) = (iat) 


Solving for 2 and 2!: 


7) i (W) 
a 
+ M2 
ee oe 

>. k fy (W) 
gt = 


m2 
£.(W) + Fi 
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Points on the curve for the W-scale are obtained by taking dif- 
ferent values of W and determining the corresponding values of z and a', 


These points are plotted anda curve may be drawn through them. 


a ; 
B adit rot va 


[J herr oxaveb bas W tc ectilav @ 


Avo & Exe bedtolaem eaintog «4 
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B. Gas-Oil Ratio Squation 
The form of the equation is: 


f, CID ak Gry, f, (W) = f), CN) (27) 


This equation is similar in form to that discussed under "A Prediction 
Equation," and is subject to the same proof, except that the scale for fy (v) 


runs in the opposite direction from the f, (W) scale: 


After substitution for x and y, this reduces to: 


ifs (‘W) = ne 7) ee —: on)- k 
1 my k-z! 2 my (x-z") 


= 


ho 


and z and z' are as before: 


mf), (W) 
f. (W) +e 
z im 


k f, (W) 
gi = i oe 


Points on the W-scale are determined by taking different values 


of W and computing the corresponding values of z and z2',. 


2 


Tet ABC be a base line, perpendicular to the three scale axes, and 
let a and b be the distances between scales. 

Draw any index line in general, FH, forming similar triangles DEF 
and DGH, where DE and GH are parallel to ABC. 


From these similar triangles: 


which rearranges to; 


© [rs 
+ 
Oo 
wo 
oO 
a 
o’ 


Substituting particular functions for x, ys, and 2: 


By a (eo 


de ee ee SEN ees 
a ar b ™ ab/(atb) 


To reduce this to the original form, let: 


a = m, 
b = Mo 
m 
ab = m3, 0r m = ee ea: 
a+b Tl + M 


This also defines the scale moduli in terms of each other, showing 
that after m, and M, are selected, ue) may be computed from them. 

To find the position of the center scale in terms of distance between 
the outside scales, note that the distance from A to B in terms of AC is; 


a 


a+b (AC) 


or 


Sie oie) 
chee 
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D. Average Gas-Oil Ratio Equation 


The equation form is: 


or 


fy (u) + fy (v) = £, @™) 


The proof is identical to that for the liquid saturation equation 


nomograph. 


hh 


APPENDIX V 
Data Concerning Reservoir "X" 


A Mid-Continent Limestone Reservoir 


Basic Data 


Average Depth 5300 =f i,. 
Original Reservoir Pressure iis pete 
Reservoir Temperature 120° FF. 
Average Connate Water Saturation 33.4% 

e Original Formation Volume Factor, & eo 1.280 
Original Solution Gas-Oil Ratio 500 = 


P.V.T. Data 


Pressure Ss ’ 

(psig) (SCF/Bb1 .STD) ey (Ours) 
1500 500 1.280 ees 
1250 118 1.260 32 
1000 392 1Z2h0 ig arad 

750 336 1216 iol 
S500 278 1.193 Hee 
250 208 ao 2n0 

100 145 1.136 Za3t 
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Volumetrically Weighted Average Bottom Hole Pressures 
Cumulative Oil 


BHP Produc tion 

Date (Psig) (Bbis.) 
Original Lii3 ate 
10~3-52 165) 152,330 
4-13-53 1558 3535995 
Tees 1517 13,999 
10-1-53 1473 5435783 
1--5h Wh 631, 76h 
7-6-5) 132 866, 063 
1-3-55 1209 1,086, 797 


Reservoir Fluid Characteristics: 


Gravity of stock tank oil h1s5° APT at 60° F, 
Saturation pressure 1487 PSIG at 128° F, 
Solution gas-oil ratio 48h Std. Cu.Ft./Resid. Bbl. 
Formation volume factor 1.28 Res. Bbl./Resid. Bbl. 
Viscosity of reservoir oil 1.28 Centipoise at 1487 PSI + 128°F, 
Analysis of Reservoir Fluid Sample 

Weight % Mol % 
Methane and Lighter ao. 26.35 
Ethane 1.70 706 
Propane het S 10.62 
Tso-butane 0.82 . ett 
N-butane 2555 5 46 
Iso—pentane Pest 2.38 
N-pentane 1,52 2 63 
Hexanes 2.81 4, .07 
Heavier 82 .08 39.66 


100.00 100.00 


18 
Gas Viscosity 


Pressure 


(Psig) 
1300 
1100 

900 
700 
500 
300 


100 


Gas gravity at 60°F. + 14.7 psi = 0.819 (Air 


0.01477 
0.01401 
0.01355 
0.01280 
0.01228 
0.01193 


0.01158 


Reservoir Gas Analysis 


Analysis of Gas Sample 


Methane 
Ethane 
Propane 
N-Butane 
Iso-Butane 
N-Pentane 
Iso-—Pentane 
Cyclo-Pentane 
Hexanes plus 
Nitrogen 
Oxygen 
Helium 
Hydrogen 

CO5 


Ho5 


Oe. 
OoL 


Orel 
100 .00 


1.000) 


hé 


Relative Permeability Data 


7 


Sy, Ky 
(4) Ke 
98 0.01 
eal O.OLSL 
96 0.033 
95 0.057 
90 0.373 
85 1,16 
80 3.20 
(6 8.0 
Summary of Production: 
Producing Production 
Wells at End Oil (Bbis. Gas (MCF 
Year of Year Year Cumulative Year Cumulative 
1951 5 274295 27,295 13,655 13,655 
1952 26 1855930 2135225 93,105 106, 760 
1953 29 414,034 627,259 254,261 361,021 
195k 31 456,238 1,083,497 753,605 1,114,626 


Production prior to bubble point: 575,000 bbls. STO 


Total Volume of Original Oil in Place: 
N = 18.000,000 Bbls. STO (Calculated by Material Balance and 
verified by volumetric value) 
On Basis of N = 1 Bbl. STO, 
Oil Production, B.P. to 1209 psig 
n = 0.0319 Bbl. 
Gas Prodiction, BP to 1209 psig 
G = 46.75 SCF 
Instantaneous Gas-Oil Ratio at 1209 psig 


rj = 2400 SCF/Bbl. 
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K, = Sy = 0.33 
Scale Labeled 
Label Limits 
P 1300 to 
100 
12" 
apart 
n O.03 
O50 
6" 
from others Sy 
Limiting Values of 5;: 
P = M500 s ny = ©2035 
P = LOO; ny = orale) 
Using 
(1-n, ) 
= x a 
s, S + (1 Py eae 


Oo 


Ko 
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J E Oo 
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1000 
900 
800 
700 
600 
500 
00 
300 
200 
100 


ee 

1-Siy 

Bo 
0.5203 
0.5203 


Limits 


1.264, 
1.135 
0.129 
l=nj 


“One! 


0.90 


Computation of Scale Lengths for Liquid Saturation Equation 


Log (S, - K,) = Log (K, 6 ) + Log (1-n;) 


= 0.5203 


Function 


log Ky, 


log(1-n;) 


log(S,- 


(3) 
Ky & 


0.6577 
0.6535 
0.6493 
0.6452 
0.610 


0.6358 ; 


O,o311. 
0.6264, 
0.6207 
0.6155 
0.6093 
0.6015 
0.5905 


(3) 
(1)x(2) 


0 06379 
0.5315 


K,) 


Range of 
Function 


9 .81803-10 
9.77122-10 
0.04681 | 


9.98677-10 
9 95h2h-10 
0.03253 


(i) 
Log K,8 


9 .81803-10 
9.81525 
9.812h5 
9 «80969 
9.80686 
9.80332 
9.80010 
9.79685 
9.79288 
9.18923 
9.7883 
9.7792 
$7 7122-10 


(4) 
(3) +Sy7=S, 


0.9719 
0.8655 


Actual 
Scale Scale 
Length Modulus 
1.03" 300 
9.7159" 300 
150 
(5) (6) 
300XlogKko8 (5)-» 
295 09 14.03 
29.575 13.209 
29443 «735 12 369 
29,2 .907 11.541 
292 .058 10.692 
291,0.996 9.630 
291,0.030 8 66h, 
29392055 7.689 
2937.86) - 6.498 
2936.769 5.03 
2935.49 083 
2933.772 2 1,06 
2931.366% 0 
Limit Values of 
57K, LogS;-K, 
0.6379  9.8080-10 
0.5315 9.725146-10 


Plotting 
Modulus 
300 
300 
150 
(iE) (2) 

n (1-n) 
0.03 0.970 
0.035 0.965 
0.040 0.960 
O.0l5 0.955 
0.050 0.950 
0.060 0.90 
0.065 0.935 
0.070 0.930 
0.075 0.925 
0.080 0.920 
0.085 0.915 
0.090 0.910 
0.95 0.905 
0.20 0.90 

Log x 150 
170.720 
158.819 

11,901 

Scale 
Length 


scale 


Equation 


x=300 log K,, B 


y=300 log(1-n;) 


). 
log(1l=n) 


9.98677-10 


9.98153 
9.98227 
9.98000 
9.97772 
9.975h3 
9.97313 
9.97081 


9.968 48 


9.96614 
9.96379 
9.9612 
9.9590) 
9.95665 
9.952) 


Actual Log 


Sp-K 


-~0.19520 
~0.27454 


(41) 
300X(3) 


2996.031 
2995 .359 
299.681 
299, .000 
2993.316 
2992 .629 
2991 6939 
2991.2h3 
2990 .5hh 
2989.82 
2989 .137 
2988 6426 
2987.712 
2986.995 
2986.272 


Actual Log 
x150 
-29.280 
-41.181 
-11.901 
Diff. = 
scale 
Length 


(5) 
(L) = 


9.759 
9.087 
8.409 
{3126 
7 Ob 
6.357 
5.667 
4.971 
272 
3.50 
2 865 
2.15) 
1.Uh0 
Ocites 
®) 
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